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Abstract. Recent X-ray observations have established 
that collisions between subclusters of galaxies are rather 
common phenomena. Prompted by such observations, we 
have performed N-body simulations of two equal-mass 
subclusters of galaxies, which are going to merge. We first 
have confirmed that only a part of kinetic energy associ- 
ated with the relative motion of two subclusters is con- 
verted to internal energy of each subcluster at the first 
encounter. Two subclusters, therefore, once pass through 
each other, turn around due to gravitational attraction, 
and finally merge at the second or third encounter. We 
performed simulations for a variety of dark-matter (DM) 
distributions, and find that the time scale for washing 
out the double peak structures after the first encounter 
strongly depends on the distribution of the density and 
velocity distributions of the DM. It takes longer when the 
DM is spatially extended and/or when the velocity dis- 
tribution of the DM has a Gaussian shape, rather than a 
uniform distribution. According to our calculation it takes 
more than 4 x 10 9 yr after the first encounter until the den- 
sity contour shows only a single peak. 

To explain the high fraction of the clusters with 
substructures among nearby clusters, Richstone,Loeb & 
Turner (1992; hereafter RLT) required recent(~ 10 9 /i~V) 
cluster formation. However, our results show that the 
timescale for subcluster merging is still uncertain and pos- 
sibly much longer than the time scale assumed by RLT. 
Caution should be exercised when concluding that the 
density of the universe is high by using RLT's method. 

Key words: Galaxies: clustering - dark matter - large- 
scale structure of Universe 
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1. Introduction 

The existence of double or multiple peaks in the spatial 
mass distribution has been established for a number of 
clusters of galaxies. Jones & Forman (1990), for example, 
have reported that about 30 % of 250 clusters of galaxies 
have double or multiple peaks in the X-ray surface bright- 
ness. Since the X-ray emissivity is more sensitive to the 
density rather than to the temperature, their observations 
strongly indicate the presence of multiple peaks in the dis- 
tribution of the intracluster medium (ICM) . ROSAT ob- 
servations of the clusters of galaxies have established that 
the high fraction of the galaxy clusters have double or 
multiple peaks in the distribution of the X-ray emitting 
gas (e.g. Briel et al. 1991). The existence of such double or 
multiple peaks has been also known for the galaxy distri- 
bution in ~ 40 % of 40 clusters of galaxies (Geller & Beers 
1982). Moreover, Miralda-Escude (1991) has recently in- 
ferred from the analysis of the gravitational lensing clus- 
ters that the distribution of the gravitational mass, that is 
the mass of the all constituents including dark matter, in 
the clusters of galaxies is not spherically symmetric, but 
has multiple peaks. These observations of the distribution 
of ICM, galaxies and dark matter in clusters suggest fre- 
quent occurrence of merger events between subclusters in 
clusters of galaxies. 

Frequent collisions between subclusters are of great im- 
portance for measuring the geometry of the universe. Re- 
garding the presence of substructuring as a signature of 
recent cluster formation, RLT estimated the lower limit to 
the cosmological density parameter to be flo ii 0.5. This 
is a rather large value and is in contradiction with the re- 
sults deduced from the another methods (e.g. Briel et al. 
1991; Yoshii & Takahara 1989). For example, the local es- 
timates of the cosmological density parameter, the baryon 
fraction in individual objects (e.g. Briel et al. 1992; White 
et al. 1993) and the results from the peculiar velocity field 
(e.g. Tully 1987), and deep galaxy number counts (Yoshii 
& Takahara 1989) support a low density universe. We thus 
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need to re-examine the basic assumptions made by RLT 
and check the significance of their obtained result. 

Their fundamental assumptions are: (1) the initial per- 
turbations are comoving with the Hubble flow, (2) the 
shape of the initial perturbations is spherically symmetric, 
and (3) the relaxation timescale, r ro i ax , on which multiple 
structures within a cluster is erased after the gravitational 
collapse of the cluster, is at most a cluster crossing time, 
of the order of ~ 10 9 /i _1 yr. The first and second assump- 
tions have been checked by Bartelmann et al. (1993). They 
introduced initial velocity fluctuation and non-spherical 
perturbation to find that RLT's conclusion gets firmer by 
these additional and more realistic conditions. 

However, the third assumption seems to be more cru- 
cial for RLT's argument since the lower limit of the cos- 
mological density parameter obtained by RLT's method 
is rather sensitive to the assumed relaxation time will be 
shown in Sec. 4. McGlynn & Fabian (1984, hereafter MF) 
have shown that it takes more than several crossing times 
after the first encounter of the substructures to erase mul- 
tiple structures after the first encounter of the substruc- 
tures unless a strong two body relaxation works during 
the collision. 

In the present study, we perform extensive studies of 
subcluster merging, following the line of MF, but adopting 
more general conditions for the initial spatial and veloc- 
ity distributions for the DM in clusters. In fact, the re- 
laxation time depends strongly on these distributions and 
the spatial distribution of the DM is still largely uncertain 
(see Sec. 2. 2). Using the simulation results, we investigate 
how the constraint on the cosmological density parameter 
using RLT's method is relaxed. In Sect. 2 we briefly de- 
scribe our model. The results of N-body simulations are 
presented in Sect. 3. The final section is devoted to sum- 
mary and discussion. 



2. Models 



Table 1. Standard model parameters. N represents the total 
number of galaxy or DM particles, M is the total mass of one 
subcluster, and m is a mass of one galaxy or DM particle. 
Indices DM and gal represent the dark matter and the galaxy 
component, respectively. 



Particles (N) 


1024 x 2 


Total mass (M) 


(8 x 10 14 M Q ) x 2 


Nv M /N gal 


4 


m D M/m ga i 


100 


M DM /M gal 


400 



2.2. Initial mass distributions in each cluster 

It is observationally known that the galaxy distribution in 
compact clusters is well approximated by the isothermal 
King distribution (King 1972), 



Pgal 



Pgal 



(rAgal) 1 



3/2 ' 



(1) 



where p® al is the central density and r ga i is the core radius 
for the galaxies in the cluster (Sarazin 1988). According 
to Bahcall (1975), we fix the core radius to r ga \ = 250 kpc 
in the present studies. The central density p° al is defined 
by J Att p ga \{r)r 2 dr = M ga i, where M ga \ is the total mass 
of the galaxies in the cluster (Table 1). Since this integral 
diverges M rmga i ~ lnr at large r, we truncate the distri- 
bution at r max = 5r ga i, which is 1.25 Mpc for r ga i = 250 
kpc. 

Unlike the galaxy distribution, the DM distribution 
is still under discussion (Hughes 1989; Fitchett 1990). 
We hence adopt the following parameterization in anal- 
ogy with the gas distribution in the cluster (Cavaliere & 
Fusco-Femiano 1976); 



2.1. Basic assumptions 

We investigate the merging of two virialized subclusters as 
an example for substructure within an unrclaxed cluster 
of galaxies. We use the same model parameters as MF to 
enable a direct comparison (see Table 1). Two subclusters 
with equal total masses are assumed, since the size of sub- 
structure is comparable to that of the main structure in 
about half of the observed clusters (Dressier & Shectman 
1988). Particle number (N) of each cluster is 1024, 20 % 
of which are assumed to be galaxies and the remaining 80 
% to DM. The mass of a dark- matter particle, muM, is 
100 times larger than the mass of a galaxy particle, m ga i. 

Here we assume that the gravitational force of the ICM 
is negligible for the evolution of the cluster. 



Pdm 



1 + {r/rv M ) : 



3a/ 2 



(2) 



where the core radius, rpM, and the slope, a, of the DM 
distribution are free parameters. The adopted parameters 
are summarized in Table 2. A constant pp M is determined 
so that Mdm = / 4:^PBM(T)r 2 dr. The combined X-ray and 
optical study of the mass distribution in the Coma clus- 
ter (Hughes 1989) has shown that the mass-follows-light 
models are the most precise to date, i.e. tdm = fgai and 
a = 1 in our notation. The study of the distribution of the 
dark matter in the A2256 (Henry et al. 1993) has shown, 
on the other hand, that the mass-follows-gas models, i.e. 
flatter distribution with a ~ 2/3, is also consistent with 
the present X-ray observations. A steeper mass distribu- 
tion than the galaxy distribution, i.e. a ~ 4/3, has not 
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Table 2. Adopted parameters for the structure of a subclus- 
ter and the resultant rebounce times (1 Gyr = 10 9 yr). The 
first column is the model names. p(r) is the density distribu- 
tion of the DM. MK, IK and PL represent the modified King, 
the isothermal King and the power-law density distributions, 
respectively, a is the slope of the density profile, fdm/ 'r rmg ai 
is the core radius of the DM distribution in unit of the core 
radius of the galaxy distribution. f(v) is the velocity profile 
of the DM. G and U represent the Gaussian distribution and 
uniform distribution, respectively. The initial relative velocity 
is given by Eq. (5), except for Model J, in which Vim = 0. 
The core radius for galaxies, r ga i, is taken to be 250 kpc, while 
-Rini = 3.125 Mpc in all models. 



Model 


p(r) 


a 


fDM/fgal 


/(«) 


Trcb(Gyr) 


A 


MK 


2/3 


0.1 


G 


9.0 


B 


MK 


2/3 


1.0 


G 


8.0 


C 


MK 


2/3 


3.0 


G 


9.5 


D 


IK 


1 


0.1 


G 


7.0 


E 


IK 


1 


1.0 


G 


8.5 


F 


IK 


1 


3.0 


G 


10.5 


G 


MK 


4/3 


0.1 


G 


6.0 


H 


MK 


4/3 


1.0 


G 


7.5 


I 


MK 


4/3 


3.0 


G 


9.0 


J 


MK 


4/3 


0.1 


G 


4.0 


Mo 


PL 


(2/3) 


(0.0) 


U 


7.0 


Mi 


PL 


(2/3) 


(0.0) 


G 


11.0 



been rejected either (Hughes 1989; Ponman 1992). Fur- 
ther more, much smaller core radii for the DM distribu- 
tion in rich cluster than the core radii of the luminous 
components are inferred from the mass distribution mea- 
surements using the gravitational lensing techniques; the 
high rate of finding arc-like structures (Le Favre et al. 
1994; Wu & Hammer 1993) and the geometrical study of 
the giant luminous arcs in distant clusters (e.g. Hammer 
1991; Miralda-Escude & Babul 1994; Wu 1994). For com- 
parison, wc also performed a simulation with a large core 
radius, i.e. 3 x r ga i. Note that MF adopted a power-law 
density profile 

Pdm = Pdm (rAgai) > (3) 

corresponding to a profile with tdm = and a — 2/3 in 
our notation. 

We assume two different velocity distributions: One 
is a Gaussian distribution with a cut-off at the velocity 
of V = ±4cr. a is determined by the condition that the 
initial system is virializcd. The second one is a uniform 
velocity distribution, from V = to V = V max , where 
Knax is determined by the condition that the initial system 
is virialized. The latter one was used by MF. The adopted 
parameters of all the calculated models are summarized 
in Table 2. 



2.3. Time evolutionary calculations of subclusters 

In a flat universe, the maximum expansion radius of per- 
turbations which collapse at a redshift z c in a f2g = 1 
universe is expressed as 

£ max = G.Oikfifa + c )- 1 /i- 2 / 3 Mpc, (4) 

where Mi 5 is the mass of the perturbation in units of 
1O 15 M . For z c ~1.0 (~ 4.0 x lO 9 /*" 1 yr from now in a 
critical density universe and even more in a low density 
universes), we find ^ max ~ 3.0/i~ 2 / 3 Mpc. If Oo < 1 or h < 
1, the maximum expansion radius of the clusters which 
collapsed at z c = 1 is larger than the above value. The 
initial separation of two clusters is set to be i?; n i = 2.5 x 
r max (= 3.125 Mpc) in all our models. Zero initial relative 
velocity with respect to the center of mass of the system , 
Vini = (model J), corresponds to a collapse earlier than 
z c = 1. Larger V[ n i correspond to the case of the larger 
maximum expansion radius, i.e. a later collapse. 
In all our calculations, except for model J, 

[2GM 

/ ]U[ x 1/2 / P. . n. -1/2 

~ 1700(km/s)(- {^-) ■ (5) 

v 1 ; V8 x 1O 14 M / V3Mpc/ v ' 

This is derived from the condition that two clusters were 
at rest at infinity. The crossing timescale, except for model 
J, is then approximately 

= -^ini/^ini 

The gravitational force is calculated by direct summa- 
tion. We use softening parameter of e — 0.6 Mpc. It has 
been shown (see p36-38 in Sarazin 1988) that the dark 
matter cannot be bound to individual galaxies but rather 
must form a continuum to prevent clusters from over re- 
laxation due to dynamical friction. In other words, the 
two-body relaxation timescale in clusters is much longer 
than the timescale of a merger. Therefore the merging of 
clusters is treated as a coUisionless process and we took a 
large value of e = 0.6Mpc which saves a lot of computating 
time. 

The Equations of motion of particles arc integrated 
by the 4-th order Runge-Kutta method. To examine the 
stability of the initial cluster configurations, we calculate 
the evolution of a single cluster for several crossing times 
(Eq. (6)), confirming no appreciable structural changes. 
The energy is conserved within 3% of the total energy in 
the every calculation. 

3. Simulations 
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3.1. General behavior 

We first calculate models with a = 2/3, the same density 
profile as those of MF, pdm cx r~ 2 at r ^> r ga i. The core 
radius is assumed to be rDM/^gai = 0.1, 1.0, and 3.0, for 
models A, B, and C, respectively (see Table 2), and the 
initial velocity is given by Eq. (5). Figure lEl shows a three- 
dimensional distribution of the DM particles for model B. 
Two clusters pass through each other, turn around due to 
mutual attraction, and recollide. 



motion. The half-mass radius increases right after the en- 
counter. This indicates that the energy associated with 
the relative motion of the clusters is partly transferred 
into internal energy within the clusters during the first 
encounter. Nevertheless, further encounters are necessary 
to damp the oscillating motions of two subclusters com- 
pletely, by forming a big cluster. 



Fig. 1. Three-dimensional distribution of the DM particle po- 
sitions for model B. The size of the box is 3 Mpc x 3 Mpc 
x 8 Mpc. The elapsed times in unit of t rm cr are, from top 
to bottom, t — 0.0,1.6,3.2, and 4.8 in the left panel, and 
t = 6.4, 8.0, 9.6, and 11.2 in the right panel. 

The upper panel of Fig. 2 illustrates the time vari- 
ations of the distance of the center of mass of the two 
clusters in units of the crossing time, t a (Eq. (6)). The 
lower panel of Fig. 2 displays the half-mass radius of the 
DM in a cluster evolving with time, which is a rough in- 
dicator of the kinetic energy associated with the internal 

1 The number of particles in Fig. 1 is reduced for the sake of 
saving the size of PS file in this online preprint. 



Fig. 2. The relative distance between the centers of the DM 
distributions of the two subclusters (the upper panel) and in 
the half-mass radius of the DM in one cluster (the lower panel) 
for models with a = 2/3 versus time. In both panels, the solid 
line corresponds to the model with tdm = 3r ga i, the dotted 
line to tdm = r ga i, and the dashed line to tdm = 0.1r ga j, 
respectively. The units of time and distance are t c (Eq. (4)) 
and 1 Mpc, respectively. 

Our main aim of this paper is to investigate the ques- 
tion how long the bimodal structure is preserved after the 
gravitational collapse of the cluster. The time of the first 
encounter of two substructures in our simulations corre- 
sponds to the cluster collapse. Since at least two encoun- 
ters are necessary to wash out the bimodal structure in our 
all simulated models, the time interval between the first 
and second encounter, which we define as the rebounce 
time, gives a lower limit to the merging time-scale after 
the gravitational collapse. The resultant rebounce times 
are summarized in Table 2. We see that it is insensitive to 
the core radius of the DM in the models with a = 2/3. 



F. E. Nakamura et al.: Subcluster merging in clusters 



of galaxies and the cosmological density parameter 



5 



3.2. Comparison with MF 

For comparison with the results of MF, we performed sim- 
ilar simulations in model M . The density profile is power- 
law (Eq. (3)) (corresponding to the models with a = 2/3 
and ruM = 0) and the velocity distribution of the DM is 
uniform. Figure 3 displays the distance of the two subclus- 
ters and the half-mass radius, respectively, as a function 
of time. The rebounce time of model Mo is 7 x 10 9 yr. Since 
the time unit defined in MF corresponds to 3 x 10 9 yr in 
our models, this rebounce time corresponds to 2.3 in MF's 
time units. This is similar to the rebounce time of model A 
in MF which is the most similar to our models. Therefore, 
our simulation code could well reproduce MF's results. 



Fig. 4. The same as Fig. 2 but for a = 1. 



Fig. 3. The same as Fig. 2 but for the models without cores. 
The dotted line corresponds to MF's original model with a flat 
velocity profile (model Mo), while the solid line represents the 
model with the Gaussian velocity distribution (model Mi) 



To investigate the effect of velocity distribution of the 
DM, a more realistic velocity distribution, a Gaussian dis- 
tribution, is used in Mi, that is the Gaussian distribution. 
The rebounce time is about 1.5 times longer than MF's 
model M . 



3.3. Models with a = l and 4/3 

In models with steeper density profiles, the rebounce time 
is rather sensitive to the core radius. Figure 4 is the same 
as Fig. 2 but for models with a = 1; the solid line repre- 
sents a model with tdm = 3r ga i, the dashed line a model 
with tdm = 0.1r ga i. The rebounce time responds with the 
core radius much more than the case with a = 2/3; i rc b is 
reduced by a factor of 1.5, when we decrease tdm by one 
order of magnitude. A similar trend can be seen for the 
case of a = 4/3 (see Table 2). The core-radius dependence 
becomes weaker as a increases. 

3.4- Bound models: dependence on Vi„i 

As discussed in Sect. 2.3, the case of the zero initial rel- 
ative velocity correspond to the case in which the clus- 
ter collapses at earlier than z c = 1. To demonstrate that 
multiple-peak structures still persist in the present mass 
distribution even if clusters collapsed at such an early 
epoch, we calculate a model with Vi„; = 0. Since the re- 
bounce time of model G is the shortest among our calcu- 
lated models, we set the same initial conditions as model 
G except V\ n { in model J. The rebounce time for models 
with Vini = but with other combinations of a and Tdm 
should be longer than the result of model J. 

The time variation of the separation is illustrated in 
Fig. 5. Clearly, the rebounce time decreases by a ~ 30 %, 
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however, still it is factor Ah longer than the relaxation time 
assumed by RLT, that is 10 9 /i _1 y r - Further oscillatory 
motion with an amplitude of IMpc still remains after the 
second encounter. This result shows that we do not have 
to require the recent formation of nearby clusters with 
substructures as RLT did, where they assumed the sub- 
structuring clusters have collapsed within 10 9 ft. _1 yr from 
now. 



Fig. 5. The same as upper panel of Fig. 2, but for a = 4/3 
and roM/fgai = 0.1. The dashed and solid lines represent the 
case with non-zero and zero initial velocity, respectively. The 
unit of time is 10 9 yr. 



Fig. 6. Density distribution of some calculated models. The 
upper panels shows models with a = 2/3 together with model 
MO, while the lower panels shows models with a = 4/3. In both 
panels, the solid line corresponds to the case with vom = 3r ga j , 
the dotted line roM = r ga i, the dashed line vum = 0.1r ga i, and 
the dash-dotted line MO, respectively. 



3. 5. Interpretation of the results 

Here we try to give the physical interpretations of the 
above results. As noted in Sect. 2.3, the two-body relax- 
ation is not effective for subcluster mergers. Instead, vio- 
lent relaxation, an energy exchange through variations of 
the potentials, plays a dominant role. The efficiency of the 
energy conversion then critically depends on how the po- 
tential changes with time, which, in turn, depends on the 
precise DM distribution in the cluster. The specific energy 
of a particle, e = v 2 /2 + 0, varies as 



de dv d<j) _ d<j> 
di =V di + dt +VV * = dt> 



(7) 



where v is the velocity of a particle, <j> is the gravita- 
tional potential, and we used equation of motion, dv/dt = 
—V<j> (Lyndon-Bell 1967). During the subcluster collisions, 
moreover, we may approximate as d<f>/dt ~ (dcf>/dx)v. 
That is, the steeper is the potential curve, the quicker 
proceeds the merging process due to an efficient energy 
conversion. 

In Fig. 6 we illustrate the density profiles of some 
calculated models, whereas Fig. 7 displays a potential 
of each model cluster. For the same core radius, say, 
'"DM/^gai = 0.1, the potential curve is steeper, when a 
is larger, so that the relaxation time becomes shorter for 
models with larger a values. When a cluster has a rela- 
tively large core radius, ruM, however, differences in po- 
tential curves are not significant, which explains why i re b 
does not change in models with roM/^gai = 1.0 or 3.0. 

Likewise, for the same a values, say, a = 4/3, the po- 
tential curve is steeper, when tdm is smaller. The relax- 
ation time is thus a decreasing function of the core radius 
when a = 1 or 4/3. In the models with a = 2/3, on the 
other hand, this tendency is not clear, because the poten- 
tial gradients do not differ so much from each other. 

We also anticipate that the relaxation time is shorter, 
when the fraction of high-velocity particles is large. Fu- 
nato et al. (1992) have shown that particles with initially 
higher energy get more energy than initially lower energy 
particles during the so called violent relaxation process. 
In this sense, a flat velocity distribution is favoured than 
the Gaussian distribution in terms of triggering a rapid 
relaxation. The distinction between models M and Mi is 
thus explained. 
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riod after the first encounter. Our results demonstrate that 
the rebounce time (hence the relaxation time) depends on 
the precise density and velocity distribution of the DM. 
There are two key parameters for the density profile, the 
core radius and the power index of the DM distribution 
in the outer region. We find that a smaller core radius 
or a steeper DM distribution yields a quicker relaxation. 
This is understood as the result of violent relaxation, in 
which a more efficient energy conversion from the relative 
motion to the internal motion occurs when the potential 
is steeper; namely, when a core radius is smaller and the 
density distribution is steeper. Our results show that it 
takes more than 4/i times longer than the time scale as- 
sumed by RLT to erase substructures in the clusters after 
the first encounter even in the most quick merging cases, 
if the cluster collapsed later than a redshift of 1. 



Fig. 7. Gravitational potentials generated by one cluster for 
each mass distribution model. Lines are same as Fig. 2. 



3.6. Non-central collisions 

So far only head-on collisions of subclusters were con- 
sidered. We have also performed the simulations for an 
impact parameter b = 250 kpc. We find an even longer 
rebounce time than in central collision. This can be ex- 
plained as follows: For b ^ 0, the center of each cluster 
does not pass through the center of the other component 
where the potential is the deepest. This means, each parti- 
cle experiences a smaller potential variation than the case 
with b = 0. Since the energy-conversion efficiency from 
the bulk motion to the internal motion depends critically 
on the potential changes (Eq. (7)), the relaxation times 
are systematically longer when 6^0, compared with the 
models with 6 = 0. Therefore, a head-on collision yields 
a shorter relaxation time than an offcenter collision. The 
possibility of non-zero impact parameter in real cluster 
collision (Fabian & Daines 1991) strengthens our conclu- 
sion mentioned in Sec. 3. 4. 

4. Discussion and conclusion 

We have performed N-body simulations to study the re- 
laxation process of two, equal-mass colliding subcluster of 
galaxies. We have demonstrated that it is difficult to erase 
the bimodal structures in the dark matter distributions 
at the first encounter and that it remains for a long pe- 



Fig. 8. Relation between flo and 5t for various values of SF, 
the fraction of clusters formed between To — St and To with To 
being the present age of the universe (cf. RLT) . 



This large merging time scale results in an appreciable 
decrease in the lower limit of S! . Figure 8 displays the 
relation between fio and St for various values of SF, which 
represents the fraction of clusters formed between To — St 
and T , where T is the present age of the Universe (cf. 
RLT). We find an approximate fitting formula from Fig. 



n ;> 



0.328F 



5t/T {n ,X ) 



(8) 



for 0.25 <S SF £ 0.40 and 0.0 < St/T & 0.6. Numerically 
solving Eq. (8) for Q yields an approximate formula, 



n ;> 



0.26.5T 
St H ' 



(9) 



which is almost independent on Ao. It shows that the de- 
rived lower limit of Oo is strongly depending on the as- 
sumed relaxation time. RLT took SF <; 0.25 mainly from 
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observed X-ray surface brightness contours by Jones & 
Forman(1990) and derived fi £ 0.5 for St £ O.Lff^ 1 . Our 
simulations demonstrate that the rebounce timescale is 
ireb/*RLT ~ 4 - 10ft, where £rlt = 10 9 ft _1 yr is merging 
timescale adopted by RLT. Since the relaxation time is 
longer than the rebounce timescale, we have 

St £ t rcb ~ (4 - 10) x 10 9 yr, (10) 

as the relaxation time of the DM distribution. If the gas 
always follows the DM distribution, we can apply this re- 
sult to estimate the lower limit of Qq with this substruc- 
ture frequency deduced from X-ray observations. It shows 
that even if St = i ro b, the derived lower limit of the density 
of the universe is more than a factor of 2 smaller than the 
value derived by RLT as 

O^O.lS^pi)" 1 . (11) 

If we take into account the fact that the oscillatory 
motion of merging substructures still survives after the 
second encounter and that the adopted value is obtained 
by using the model with an extremely centrally concen- 
trated mass distribution which leads to the shortest relax- 
ation time, the true relaxation time could be longer than 
this and a much lower density universe could be consistent 
with the observed substructure frequency. 

Finally, we would like to note the applicability of our 
present results to the estimation of the mean density of 
the universe based on the statistics of the appearance of 
the substructures in clusters. In the above discussion, we 
have assumed that the morphology of the X-ray emitting 
gas always follows the morphology of the dark matter dis- 
tributions. However, it is still unclear whether the hot gas 
follows the DM or not. There is evidence which supports 
this assumption both in simulations (Roettiger et al. 1993) 
and in the real world (Briel et al. 1992). Roettiger et al. 
(1993) have shown that even after the encounter of two 
clusters the X-ray emitting gas almost traces the dark 
matter distributions and they have suggested that the 
substructure in the Coma cluster found by X-ray mor- 
phology (Briel et al. 1992) is in this post encounter stage. 
However, there is also a counter evidence against this as- 
sumption in simulations (Evrard 1990; Schindler & Miiller 
1993) where the bimodal structure of the X-ray emitting 
gas disappears when two substructures meet although the 
dark matter of the two substructures passes through each 
other (Evrard 1990; Schindler & Bohringer 1993). Proba- 
bly the both cases are occurring in real clusters depending 
on their merging conditions but it is not sure which is the 
dominant process. If the fraction of merger events rep- 
resented by the former case is dominant, we can safely 
assume the above mentioned assumption. However, if the 
fraction of the merging events represented by the latter 
case is dominant, the frequency of substructures in dark 
matter distribution deduced from X-ray data is highly un- 
derestimated. Fortunately, there is now a new useful tool 



to study the dark matter distribution in clusters directly, 
that is gravitational lensing technique (e.g. Fort & Mel- 
lier 1994). If the latter is the case, we will find a large 
discrepancy of substructure frequency between the results 
deduced by the gravitational lensing technique and the X- 
ray data, that is the former will be much higher. This will 
be a good test which of above two processes is dominant 
process in real world. 

The special nature of the initial conditions in our sim- 
ulations would also limit the applicability of our results 
to the real cluster merging. There have been extensive N- 
body calculations based on somewhat more cosmologically 
motivated initial conditions than are the present calcu- 
lations (e.g. Evrard et al. 1993; Crone, Evrard & Rich- 
stone 1994; Jing et al. 1994). Although their approaches 
are basically more straightforward to compare their re- 
sults with real cluster observations, their present results 
are still largely limited by their numerical resolution, size 
of the simulation boxes, number of the simulated samples, 
baryon fraction, and so on. Actually, the results obtained 
by this approach are still controversial. Jing et al. (1994) 
have obtained much a larger fraction of clusters having 
substructures in a low density universe than Evrard et 
al. (1993) who found that the morphology of clusters in 
low-density universes are much more regular, spherically 
symmetric than those in the critical density universe. In 
this situation, our approach has a complementary mean- 
ing to examine the cluster merging process. Our present 
study is enough to make a caution to require recent cluster 
formation to explain the high fraction of nearby clusters 
having substructures and to conclude that the density of 
the universe prefers to be high by using RLT's method. 

In conclusion, the result of a high density universe RLT 
derived should be still in question because multiple X-ray 
clusters cannot be easily considered as evidence for recent 
cluster formation according to our simulation results. Fur- 
ther studies to clarify the dependence on the relaxation 
timescale is required. 
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